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1. Introduction 

For the next generation of Extremely Large Telescopes (ELTs), several Wide Field Adaptive 
Optics (WFAO) concepts like Ground-layer AO (GLAO), Multi-conjuguate AO (MCAO), 
Laser Tomography AO (LTAO) or Multi-Object AO (MOAO) have been proposed. These 
systems will play a crucial role as most of them could be used as first light instruments. 
It is then of prime importance to understand the fundamental limitations, the possible 
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optimizations, and the expected performance of those systems. The goal of this paper 
is to bring a quahtative and quantitative analysis of the most general problem of phase 
reconstruction and correction in multi-Guide Stars (GSs) AO. We try to point out some 
general trends shared by all these future systems and we propose physical interpretations of 
the results. 

All these WFAO concepts have in common that they need the knowledge of the 
3-dimensional turbulent volume. The method of measuring the 3-dimensional atmospheric 
phase perturbations was proposed by Tallon and Foy [1]: light from several Guide Stars 
(GSs) is used to probe the instantaneous 3-dimensional phase perturbations, the turbu- 
lent volume is then reconstructed by solving an inverse problem. This technique, called 
atmospheric tomography, was later improved by Johnston and Welsh [2], Ellerbroek [3] 
or Fusco [4] for instance. The tomographic reconstruction error only depends on the GSs 
constellation, atmospheric conditions and Wave Front Sensor (WFS) characteristics. The 
geometry of the system set the first two fundamental limitations of the tomography: unseen 
modes and unseen turbulence. Unseen turbulence that originates from partial pupil overlap 
can play a non negligible role in the global residual variance (e.g. [5,6]), however, in the 
frame of ELTs this term becomes negligible [7] . We then choose to neglect the error due to 
unseen turbulence in the present study. Instead, in a first part of this paper, we focus on 
unseen modes and their impact on the phase estimation depending on the reconstructor 
choice. 

In addition to this fundamental limitation, in a second part of this paper, we investigate 
an additional term due to model/statistical errors. This term comes from a mis-knowledge 
of system and atmospheric conditions, which limits the tomographic reconstruction. For 
instance, Conan et al. [8] and Tokovinin et al. [9] have studied the sensitivity of the 
reconstruction error to an error on the Cn2 profile; LeLouarn and Tallon [10] have explored 
the influence of the reconstructed turbulent layers altitudes compared to the real turbulent 
volume; Fusco et al. [11], the effect of a reconstruction on a limited number of Equivalent 
Layers (EL). All these studies suggest that the tomography is particularly robust to 
model/statistical errors. We generalize these approaches for the ELTs, and try to give a 
comprehensive study of the impact of these model errors. 

Once the 3-dimensional turbulent volume is reconstructed, different types of corrections 
can be applied. The correction step consists in a projection of the reconstructed tomographic 
phase on the Deformable Mirrors (DMs). It only depends on the number and position 
of the DMs with respect to the reconstructed layers, and the optimization area. If the 
corrected area is not increased compared to classical AO, only one DM conjugated to the 
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telescope pupil is required to perform the correction (MOAO, LTAO). The main advantage 
of those methods is that any target selected in the FoV can be corrected: the sky coverage is 
significantly increased. Moreover, there are no additional terms of errors due to projection. 
For a larger corrected field, a trade off between the optimization area and the DMs position 
is applied. This additional error term was defined by Rigaut et al. [12] and Tokovinin et 
al. [13] as generalized fitting. For instance, with one DM conjugated to the telescope pupil 
(GLAO), a wide field can be corrected uniformly, however the performance is dramatically 
limited by this error term. To reduce the impact of the generalized fitting, one must use more 
DMs optically conjugated to the turbulent volume. By correcting the turbulence volume 
above the telescope, the efficiency of the correction is significantly improved (MCAO). In 
the third part of the paper, we investigate this particular error term and its impact for 
ELTs. We compare the relative performance expected for the different WFAO concepts, we 
explore the effect of different reconstructors, and we discuss the choice of the number of DMs. 

We choose to conduct all these analyses in the Fourier spatial frequency domain because it 
allows us to perform a fast investigation of the parameter space, and it brings a comprehen- 
sive view of the impact of the different terms. The use of spatial frequency domain techniques 
was initially proposed by Rigaut et al. [14], who derived the analytical expressions for the 
five most fundamental limits of any Natural Guide Star (NGS)-based AO system: fitting 
error, angular anisoplanatism, servo-lag, WFS noise and WFS aliasing. This single conju- 
gated NGS-AO case was later extended by Jolissaint et al. [15], who included the correlation 
between anisoplanatism and servo-lag error, and generalized the formalism for 2-dimensional 
systems. In the other hand, Tokovinin et al. [9,13] have proposed an approach for evaluating, 
respectively, the wavefront fitting error with a limited number of DMs in a MCAO system 
(the so-called generalized fitting), and the tomographic reconstruction error from a limited 
number of noisy measurements (see also [16]). Following their approach. Gavel [17] has gen- 
eralized the optimal tomographic reconstruction for spherical waves. Although this approach 
requires some simplifications, it allows to treat the so called "focus anisoplanatism" error. 
Finally, Ellerbroek [18, 19] has derived an integrated approach that takes into account the 
fundamental error sources and their correlations for the general case of multi-GSs AO and 
MCAO. 

We follow the Ellerbroek's approach, with additional developments to account for model 
errors and a different formalism. For sake of simplicity, we choose to not include temporal 
behaviors, aliasing effects or closed loop considerations. In addition, we will only consider 
NGS. This NGS open-loop hypothesis could appear to be restrictive, however, it allows us 
to obtain simple analytical formulas. Considering the complexity of WFAO systems, it is 
instructive to isolate specific factors in order to disentangle their individual effects. All these 
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terms can later be included in a second step, or for a real system design for instance. 

In section 2 we briefly summarize the theoretical basis of Fourier simulation and we apply 
this method to WFAO In Sect. 3. Sect. 4 describes the simulation cases used thorough 
the paper. Sect. 5 is dedicated to the pure tomographic reconstruction and the impact of 
unseen frequencies. Sect. 6, investigates the impact of model errors on the tomographic 
phase reconstruction and Sect. 7 the impact of statistical errors. Finally, Sect. 8 focuses on 
the phase control and projection on DMs with some examples of different WFAO systems 
and their expected performance for a 42m telescope. Conclusions are given in Sect. 9. 

2. Fourier AO Modeling 

2. A. Introduction to Fourier simulations 

The starting point of the Fourier method is to assume that the optical system (phase 
propagation, WFS measurements, DMs commands) is linear and spatially shift-invariant. In 
that case, all the usual operators are diagonals with respect to the spatial frequencies and 
simply act as spatial filters in the Fourier domain. Therefore, each equation can be written 
frequency by frequency. For instance, the tomographic phase reconstruction is derived and 
evaluated one Fourier component at a time. In addition, our purpose is to derive regularized 
reconstructors which make use of the knowledge of the phase and noise statistics. Therefore, 
these statistics must also be described by linear shift invariant spatial filters. This is 
achieved by assuming gaussian and stationary statistics. In that case, the phase and noise 
statistics can be entirely characterized by their 2nd order moments, more precisely, by 
their Power Spectrum Density (PSD). As a consequence, the residual phase in the pupil 
for one direction of interest is itself fully characterized by its PSD. We will make use of 
this residual phase PSD to carry our analysis. This is similar to the usual characterization 
of the turbulent phase by Kolmogorov / Von Karman PSD, except that we derive a PSD 
for the residual phase. Similarly to what is usually done for the turbulent phase, we can 
then construct from this residual phase PSD (i) the phase structure function, the optical 
transfer function, the long exposure PSF, (ii) compute the residual phase variance, (iii) 
generate instantaneous corrected phase screens and instantaneous AO corrected PSFs. The 
derivation of the residual phase PSD for WFAO system is described in Sect. 3. 
The main advantage of the Fourier method is that the computational complexity is 
significantly reduced compared to spatial domain modeling. The main limitation is that 
aperture-edge effects and boundary conditions, which cannot be represented by shift- 
invariant spatial filters, are neglected. Therefore, the Fourier modeling only applies to the 
ideal case of infinite aperture systems. Note that this only applies for the calculation of the 
PSD, as for performance evaluation we will be able to account for a finite aperture (e.g. 
2.B). The main assumption is then that all the effects of incomplete beam overlap in the 
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upper atmospheric layers are neglected. The impact of such an assumption is discussed 
in [7], based on a comparison with a full E2E simulation code. Briefly, they show that as 
long as the different beams are superimposed in the upper layer, the error due to unseen 
turbulence has no impact on the tomographic reconstruction. For a 42m telescope and a 
typical atmospheric profile, this would correspond to a maximal GS diameter constellation 
around 8 arcmin. 



2.B. From PSD to Variance and Strehl Ratio 

In the following, we will characterize the performance in term of residual variance and Strehl 
Ratio (SR). Following EUerbroek [19] and Jolissaint [15], we define the residual phase variance 
as the piston-filtered integral of the residual phase PSD. This reads: 

oo 

aL,, = / FSBrimim (i) 



with: 



2Jii7rDfy' 



(2) 



nDf 

where D is the telescope pupil diameter. In all this paper, D is set to 42m. 

The SR is derived from the long-exposure PSF, which is computed from the PSD as described 

in Jolissaint et al. [15] 

3. Fourier WFAO modeling 

The goal of this section is to derive the residual PSD when dealing with WFAO systems. 
3. A. System description and notation 

The atmospheric profile will be divided into Nl discrete independent layers located at alti- 
tudes {hn}- Each turbulent layer is described by its own power spectrum C^^, that can be 
written for the n^^ layer as: 

C^„ = A„C*0 (3) 
is a Von-Karman power spectrum defined by: 

/-[ \ 5/3 . 1 ^ -11/6 

C,(f) = 0^023 (-) (f^ + -) (4) 

where tq is the Fried parameter and Lq is the outer scale of turbulence. In the following, we 
implicitly assume that all the layers have the same value of Lq. A„ is the fraction of turbulent 
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energy located in the n layer and is defined by: 

Nl 

\n = Cl{K)5hJY.Cl{K)5K (5) 

n=l 

where 6hn represents the width of each turbulent layer. 

Nl Nl 

We can then write that = A„,C0 = Cy,„. 

n=l n=l 

The measurement is performed with several WFSs, each WFS looking at one Guide 
Star. The number of GSs/WFSs is denoted Ngg, and the associated GS positions are given 
by a = {oi}. 

The correction is performed by Ndm DMs, optically conjugated to altitudes {h^'^^}- The 
FoV of interest, where the correction has to be optimized, is discretized into Nang angles 
f3 = {Pj}. Finally, performance is computed for NfUr direction of interests, at angles 
9 = {6k}- Figure 1 summarizes the geometry of the system. 

a| i"^ a| Tngies 



altitudes 




Fig. 1. Illustration of a system geometry. In this example, the atmosphere is 
simulated with three layers at altitudes /in=i,2,3; weighted by a Xn=i,2,3 factor 

Nl 

( A„ = 1). Two guide stars are considered in directions aj=i 2, the optimiza- 

n=l 

tion is done in three directions /3j=i,2,3, the correction is perfomed with two 
DMs at altitudes /i^if 2 and the final performance is evaluated in one direction 

e 



3.B. Residual phase, measurement equation and DM command 

We start with the calculation of the residual phase, expressed in the Fourier Domain. In 
the direction of interest 6, the residual phase ^^'^^(f) is given by the difference between the 
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atmospheric phase corrugation 06) (f) and the correction phase 0g°''^(f): 

0r(f) = 0^(f)-0r''(f) (6) 



where the tilde denotes the Fourier Transform. Assuming that propagation effects can be 
neglected, the turbulent phase resulting in the pupil when looking in the direction 9 can be 
expressed as the sum of all phase perturbations. With the Fourier shift property, it leads to: 

0,(f) = X: ^n(f)e^^-^''"^ (7) 

n=l 

where is the phase corresponding to the n^^ layer. 



The correction phase is the sum of the phase perturbations introduced by the DMs. 
This resulting correction phase in direction 9 then reads: 

Ndm 
n=l 

where ip^'^ is the correction phase done by the n*^ DM. Note that each DM is assumed to 
correct spatial frequencies up to a cut-off frequency f^^'^: 

tc - ^ (9) 

where Nact is the actuator pitch as seen from the pupil plane. Beyond this cut-off frequency, 
no correction is applied and the corresponding phase (^^*^(|f | > f^^^)=0. 



In, order to simplify the equations, in the following we adopt the notation x which refers 
to at a given frequency. For each frequency, Eq. 7 and Eq. 8 can then be written in a 
vectorial form: 

4>e = Pe <^turb 

(10) 

^corr ^ pDM -^^ 

Where cpg and (f)g°'^^ are scalar. and are line- vectors of sizes A''^,, respectively Nj^mf. 

<^turb is a column- vector of size A''^ which concatenate the Nl turbulent phases {ipn}- 
Similarly, <^dm is a column- vector of size Ndm gathering the Ndm correction phases {<^^^^}. 

Pg projects the N^ turbulent layers seen in the 9 direction in the pupil plane, and 
projects the DM phases seen in the 9 direction in the pupil. The elements of the 
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vectors Pg and Pg are simply the complex exponential appearing respectively in Eq. 7 
and Eq. 8. 

In the following we will assume that each measurement from the WFS is divided in 
two sub-measures, called hereafter measurement along x and y. For instance, this corre- 
sponds to X and y slopes when considering a Shack-Hartmann (SH) sensor. In that case, 
the most general form for the measurement is to write it as (for a given frequency): 

0— = MP^(^turb + b (11) 

Where (p"^^^ is a (2Ngs) column-vector that concatenates the measurements coming from 
the Ngs GSs. 

P^ is a matrix of size N^xNgg, where we adopt a convention "columns by rows". It projects 
the phases in altitudes coming from the Ngs ai directions in the pupil. Similarly to P^ , each 
line of the P^ matrix is filled with phase-shift complex exponential, 9 being replaced by the 
Ui direction of the considered GS. 

M is a matrix which models the WF sensing operation (e.g. Rigaut [14], Jolissaint [15], for 
an expression of the SH in the Fourier domain and Verinaud [20] for the pyramid one). It 
is a Ngsx{2Ngs) bloc diagonal matrix. Each bloc is a 2-elements column-vector modeling 
one WFS. The upper element is the WFS operator along x and the lower element is the 
WFS operator along y. Each WFS measures frequencies up to a cut-off frequency f^^^ ■ 
Beyond this cut-off frequency, the WFS operator is equal to 0. The MP^ matrix describes 
the sensitivity of the system and it will be referred as the global system interaction matrix. 
Finally, b is a column- vector of noise associated to each GS (dimension = 2Ngs). It originates 
from photon and detector noises. This WFS measurement noise is statistically independent 
for distinct GS and x and y components. The variance associated to this noise is assumed 
to be identical for all frequencies associated to a given GS, with a value denoted as af (rd^) 
for GS number i. 

Finally, one still has to express (^dm- With linear assumptions, the most general form for 
the correction phases is a sum of suitably filtered measurements defined as: 

(^DM = W0--^ (12) 

Where W is the phase volume reconstructor matrix of size {2Ngs) x N^m- It converts the 
{2Ngs) measurements into Njjm correction phases. Fusco et al. [4] have shown that this 
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matrix can be split in two independent matrices, written as: 

W = PoptWtomo (13) 

The reconstruction can thus be seen as two consecutive steps. The first one, corresponding 
to Wtomo (dimension = 2Ngs x A''^), provides an estimate of the turbulent phases on each 
reconstructed turbulent layer, it is the tomographic reconstruction. It only depends on the 
GSs configuration and the atmospheric conditions. The second step, corresponding to 
Popt (dimension = A''^, x Ndm), consists in a projection of the tomographic solution onto the 
DMs. It is a geometrical operation which provides the correction phases on the DMs from 
the tomographic estimated phase in the volume. It depends on the size of the scientific FoV, 
and on the number and positions of the DMs with respect to the reconstructed layers. In 
the following, we choose to investigate independently these two steps. 

3. C. Phase reconstruction: W-^^qj^q 

The goal of the tomographic reconstructor is to find the best possible estimate of the phase 
volume from noisy GSs measurements. The phase volume estimate is written as: 

^turb = Wtomo*^'"^^ (14) 

where the hat stands for the volume estimates. The number of sources is limited, while 
the number of turbulent layers is large not to say infinite, meaning that the problem is 
ill-conditionned and generally not directly invertible. Several methods have then been 
proposed to reconstruct the turbulent volume. We choose to investigate three reconstructors 
that are frequently used and discussed: the Minimum Mean Square Error (MMSE), the 
simple Least Square Estimator (LSE) and the Truncated LSE (TLSE) approach. 

3.C.I. MMSE reconstruction 

The MMSE estimator minimizes the residual phase variance in each reconstructed layer 
following a quadratic criterion between actual and estimated phases: 

= (ll<^turb - Wtomo(MP^(^turb + b) | ^ (15) 

In the Fourier domain, the derivation of Eq. 15 with respect to Wtomo is done frequency 
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by frequency. This leads to the optimal solution defined by: 



W 



MMSE 



-1 



fMP 



L\T, 



that can also be written as [21] 



(16) 



MP^C^„(MPJ:)T + Cb/ 

This reconstructor includes prior knowledge on the phase statistics and noise power spec- 
trum by means of C^^ and Cb- 

As the turbulent layers are independent, for each frequency, C<^^ is a diagonal matrix of size 
NlxNl. The diagonal components are given by the fraction of turbulent energy associated 
to each layer: C<^^=diag(AiC0,...,A7v^C</,). 

WFS measurement noises are statistically independent for distinct GSs and orthogonal di- 
rections, meaning that Cb is a diagonal matrix of sizes 2Ngsx2Ngs. The noise variance 
is identical for all the frequencies (white noise) with a value defined by crfi^xy) (rd^) for 
the GS number i. It can be different for different GSs, and for x and y directions: Cb oc 
diag(cT2 ^ .aly ,x ,y ) • 

If we assume the same noise for all GSs and x and y directions (Cb oc Id), we can factorize 
this noise term in the first form of w'^^^^, and re-write Eq.l6 as: 

(MPir (17) 

This form makes appear the regularization term cr^C"^ which is nothing else that the inverse 
of the Signal to Noise Ratio (SNR). We will see in Sect. 5. B the role of this regularization 
term to avoid the noise amplification due to unseen frequencies. 

3.C.2. (T)LSE reconstruction 

The LSE estimate ifturb is the one providing the best-fit to the measurements, and can be 
found by minimizing the following criterion : 

= (Il0mes - MP^Wtomo0mes|r) (18) 

The minimization of the above criteria with respect to Wtomo leads to the well known 
solution: 

(MP^)^ (19) 

It is well known and interesting to note this LSE solution corresponds to the MMSE one 
introduced in Eq. 17 for which the regularization term have been canceled (cr^C~^ tends to 



^MMSE ^ C^„(MP^)T 



w 



MMSE 



[MP^)^MP^ + a^C;„^ 



W 



LSE 



fMP!^)^MP^ 



10 



0). 

For some frequencies, the (MP^)^MP^ is not invertible or badly conditioned. For such 
frequencies, Truncated Singular Value Decomposition (TSVD) is used to truncate these 
singular values below a pre-selected threshold. For the truncated frequencies, W^^^=0 and 
(^„=0. The choice of the threshold is made empirically, by trying to find the best trade-off 
between the number of truncated frequencies and noise amplification. In simulation, it is 
relatively straightforward to find the optimal threshold with an exhaustive exploration of 
threshold levels. For each configuration, we compute the residual variance as described in 

2. B for a set of threshold levels, the optimal threshold being the one that minimizes the 
residual variance. In real systems, when accounting for overheads during on-sky operations, 
it is not easy to scan the parameter space. Hence all the following results for the TLSE case 
should be taken as optimistic ones. 

3. D. Model and statistical priors 

Few precision on the reconstructors are important at that point. 

• Firstly the number /altitudes of the reconstructed layers can be different from the 
"real" number/altitudes of atmospheric layers as introduced in Eq. 3. The choice of the 
number/altitudes of reconstructed layers comes from our knowledge of the atmospheric 
conditions. If our description of the atmospheric model is partial or wrong, we will 
commit an error. We will refer to this error as the "model error" . 

• Secondly the value of the statistical priors used in the MMSE reconstructor can also 
be different from real noise and turbulence statistics. This mis-knowledge leads to an 
additional error that will refer as the "statistical error". The TLSE does not explicitly 
include priors on phase and noise statistics. However, the truncation level is pre-defined 
based on assumed models for system/atmospheric conditions. A statistical error then 
lead to badly tuned truncation. 

We discuss the impact of the model errors in Sect. 6, and the statistical errors are inves- 
tigated in Sect. 7. 

3.E. Projection onto DMs: Popt 

The tomographic reconstruction is followed by an optimal projection of the volume estima- 
tion onto the DMs. This can be written as: 

<^DM = Popt^turb (20) 
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The filter that optimally projects the turbulent volume onto DMs depends on FoV optimiza- 
tion directions P = {Pi} and is defined by [3,4]: 



opt 



{{P^NfP^^N), - ((PrN)^P|)/3 (21) 



where, Op represents the average over all the Nang directions. N is a diagonal matrix of sizes 
Ndm^Ndm defined by N=diag(y4i,...,y4Ar^^J, with v4„=l when the considered frequency is 
smaller than the DM cut-off frequency, and otherwise. and P^. are constructed in 

the same way as P^^ and Pg in Eq. 10. 

If the mirrors positions exactly match the reconstructed layers, no projection is required 
and Popt = Id. In such a situation there is no particular optimization in the FoV, and 
no additional error due to projection. When the mirror positions do not match the 
turbulent layers, optimizing for a particular FoV position may degrade the correction in 
other directions. Trade-offs have to be made for a specific set of FoV positions, and the 
matrix Popt performs optimally these trade-offs. In such a situation, an additional error 
term due to projection arises: the so-called generalized fitting. We discuss this term in Sect. 8. 



3.F. Derivation of the residual phase power spectrum 
Putting Eq. 10, 11 and 12 in Eq. 6, we then have: 

= P^<^turb - P?^W(MP^<^turb + b) (22) 

We can now derive the residual PSD for each frequency. By definition, PSD^^* is the 
statistical average of the square modulus of the residual phase 0^*^*. Assuming that noise and 
phase are statistically independent from each other, it follows that: 

PSD^"' = 

(pL _ ppMwMP^)C^„,t(P^ - PP^WMP^)^ (23) 
+ (P?Mw)Cb,t(P?^W)T 

In the above equation, we have introduced C<^^ t? the power spectrum matrix of the phases 
in the turbulent volume and Cb,t the noise power spectrum matrix. These matrices stand for 
true noise/atmospheric statistical conditions. They are both diagonal matrices constructed 
in the same way as and Cb in Eq. 16. The differences with C<^„ and Cb, are that C^^^^t 
can be larger than C^^ if we commit a model error, and Cb,t can have different values of 
noise variance if we commit a statistical error (see 3.D). We distinguish them from the one 
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introduced as priors in Sect. 3.C.1 with the subscript t. In all the following, PSD^*^^ will be 
our working starting point. All the matrices required to construct PSD^*^* are also described 
in more details in [7]. 

4. Simulation Conditions 

Based on the residual phase PSD derived in Sect. 3.F, we want now to investigate the 
different errors introduced in Sect. 1, respectively: unseen frequencies, model/statistical 
errors and projection errors. For that we define three simulation cases representative of 
different GS/FoV configurations: 

• The first one, called "2GS" , is a simple 1-dimensional model where the atmosphere is 
only composed of 2 layers. One of the 2 layers is located in the pupil, the other one is 
at an altitude h=8km. We consider that the 2 layers have a relative turbulent strength 
profile defined by [Ai,A2]. We consider 2 GSs separated by a distance ai 2=2arcmin. 
This simple configuration is mainly used to derive didactic examples, illustrated with 
the 1-dimensional residual PSD. Fig. 2 describes the geometry of this case. 

• The second case, called "4GS", is a more realistic model where the atmosphere is 
composed of 10 layers. The altitudes and strength of the layers are summarized in 
Table 1. This turbulent volume is sensed with 4GS located on a 2arcmin diameter 
circle. All GS have the same noise variance a^=0.5 rd^. This constellation is used as a 
representative configuration for a medium FoV system. 

• The third case, called "8GS" uses the same turbulence profile, but sensed with 8GS 
located on a Sarcmin diameter circle. All GS have the same noise variance a^=0.5 rd^. 
This constellation is used as a representative configuration for a large FoV system. 



Table 1. 10 layer turbulent profile altitudes and relative layer strength. 



hn (m) 





200 


600 


1200 


3000 


4700 


7300 


8500 


9700 


11000 


An 


0.41 


0.16 


0.1 


0.09 


0.08 


0.05 


0.045 


0.035 


0.02 


0.01 



Table 2 summarizes the main parameters of these 3 simulation cases. 
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Table 2. Simulation configurations. 





2GS Config. 


4GS Config. 


8GS Config. 


Atm. profile 


21ays 


10 lays 


10 lays 


Constellation 


-X 


' — Sep=2'_' 


t 




T 


# GS 


2 


4 


8 


FoV Diameter 


2' 


2' 


5' 


Noise/GS 


0.5 rd^ 


0.5 rd^ 


0.5 rd^ 



For all simulation cases, the Fried parameter is set to ro=0.12m defined at 0.5/im and 
Lo=50m. We assume the that turbulence measurements are performed with SH WFS. In 
that case, Rigaut et al. [14] show that the measurement operation, without aliasing and 
finite exposure time, can be modeled by a M = 2j7r f^^ysmc^ndfx^y) term, where we assumed 
a square geometry for subapertures with a size (reported in the telescope pupil) equal 
to d. For our example, we consider that all WFS have the same subapertures sizes. This 
implies that all frequencies greater than fc = l/2d are not measured. We further assume 
that all the GSs have the same noise variance defined by = 0.5 rd^, where o"^ are the 
diagonal components of Cb,t- Performance is evaluated at 1.65/im (H Band) by projecting 
the volumetric residual phase in a direction of interest 6. Unless specified otherwise, it is 
evaluated at the center of the FoV {9 = 0). 

In the following sections, we make use of those three simulation cases to investigate (i) 
pure tomographic reconstruction in Sect. 5; (ii) model/statistical errors in Sect. 6 and 7; (iii) 
projection issues in Sect. 8. We want to emphasize that most of these errors strongly depend 
on the turbulence profile. More precisely, we will see that for a given isoplanatic angle, a 
relevant parameter is defined by the products between layer altitudes and GS separations: 
hn<y. In our analysis, we have chosen to fix the turbulence profile, and investigate these 
errors in term of GS separation (FoV). Numerical results presented in the following should 
then not be seen as absolute results, but rather as quantitative insights used to illustrate 
general trends. 
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Fig. 2. System geometry for the 2GS two-layers simple case. Altitude of layers 
is hi=0 and /12 = h. Two guide stars are considered in directions ai and a2, the 
separation ai2 between the GSs is 2'. Performance is evaluated at the center 
of the FoV: ^ = 0. 



5. Pure Tomographic Reconstruction 

In this section, we consider the ideal situation where noise and turbulent conditions are 
perfectly known: C^^ and Cb used in Eq. 16 are then identical to C^^^^t and Cb,t used in 
Eq. 23. This situation is generally called the "pure tomographic reconstruction". 

5. A. Analytical expressions of the residual PSD 

We first consider the 2GS case. From this simple 1-dimensional configuration, we can derive 
the analytical form of the residual phase PSD for both filters. Eq. 23 can be written when 
considering a (T)LSE as: 



PSD^l^o 



o-^ [2-cos(27r/(ai -e)h)-cos{2Trf(a2 ~e)h)] 



2sin2(7r/oi2/i)Af2 

for i < fc and for non truncated frequencies 



PSDgl^o = for f > /c or for truncated frequencies 
while with an MMSE, it reads: 

2MC^AiA20-^ [2~cos{2TTf{ai-e)h)-cos{2TTf{a2-e)h)]+a^ 



(24) 



PSD«!fQ — ^ isin^(Trfai2h)M^C'^XiX2+2M^a^C^+a-i 

[ PSD^l^o = C<^forf>/, 



for f < /e 



(25) 
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In Fig. 3 we compare these two 1-dimensional residual PSD laws for 6* = 0, as well as the 
uncorrected Von-Karman spectrum and the on-axis reconstruction for comparison. The on- 
axis reconstruction corresponds to a classical LSE Single Conjugated AO configuration and 
is computed by setting the two GSs in the direction of interest. It follows the typical law 
expected for derivative WFS [22]. Since the two GSs are on-axis, we benefit of the global flux. 



5.B. Notion of unseen frequencies and neutral frequencies 

In Fig. 3, it appears that for some specific frequencies, the tomographic residual PSDs 
(either MMSE or LSE) diverge from the on-axis noise propagation, while other frequencies 
are following the on-axis noise propagation. 




Fig. 3. Comparison of residual PSDs for the MMSE reconstruction (full line) 
and un-truncated LSE reconstruction (dashed line). Uncorrected Von-Karman 
(dotted line) and on-axis reconstruction (dashed-dotted line) are also plotted 
for comparison. In this example, both layers have the same turbulent strength 
and the WFS cut-off frequency is fc = 0.9m~^. Two unseen frequencies ap- 
pear for /i=0.215m~^ and /2=0.645m~^, two neutral frequencies appear for 
/i=0.43m-^ and /2=0.86m-^ 



To understand the behavior of these PSDs, we first consider the LSE reconstructor. From 
Eq. 24, we learn that the reconstruction error tends to infinity for each spatial frequency equal 
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to p/ {a^h), where p is an integer. Indeed, every time that the period of a frequency exactly 
match the GS separation, the sum of the phase perturbations induced by the 2 layers is 
exactly zero, the phase information for these frequencies is lost. We refer to these frequencies 
as "unseen frequencies" , as they are not sensed by the WFSs. For a given turbulence profile, 
the unseen frequencies only depend on the GS positions: the more distant the GSs, the 
smaller the spatial frequencies affected, the larger the number of unseen frequencies in a 
given frequency range of interest. The GS separation can then be understood as a spatial 
basis which sets the sensitivity of the filter to the different frequencies. 
We now investigate what is the impact of the unseen frequencies when looking in a particular 
direction of interest 6. This relies on the numerator of Eq. 24. We find that for some specific 
frequencies this numerator goes to zero. This happens each time that the sum of the phase 
perturbations induced by the 2 layers is exactly zero in the direction 6. These frequencies 
are not impacting the performance in the direction of interest, we will called them "neutral 
frequencies" . For instance when the direction of interest is at ^=0 (as in Fig. 3), these neutral 
frequencies occurs each / = 2p/{ai2h). Therefore, for this specific direction, the neutral 
frequencies are canceling one unseen frequency over two. We note that for the frequency 
exactly corresponding both to an unseen and neutral frequency, Eq.24 is undetermined. For 
that specific frequency, calculation of the limits gives: 



which is the classical on-axis propagation law [14]. 

Another example of the effect of the neutral frequencies is when the direction of interest is 
exactly lying on one of the GS (e.g. 9 = ai). In such a situation all the unseen frequencies 
are canceled by the neutral frequencies: an unseen frequency has no effect on the image 
quality in the GS directions. 

For the pure LSE reconstructor, unseen frequencies are a real issue as they lead to an over- 
amplification of noise. For these frequencies, the interaction matrix coefficients go to zero, 
direct invert is ill-conditionned and the noise is dramatically amplified. Around an unseen 
frequency, there is a set of badly-seen frequencies. For these frequencies, measurement is 
close to zero, but since the WFSs are noisy, the same problem as for the unseen frequencies 
applies. In the following, we will also refer to these frequencies as unseen frequencies. 
To avoid this noise amplification, one should use the TLSE reconstruction. In Fig. 4 we show 
the corresponding residual PSD for the TLSE reconstructor having the best threshold. Noise 
amplification is avoided, at the price of uncorrected frequencies. 




(26) 
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Fig. 4. Same as Fig.3 with TLSE reconstructor. Threshold was optimize to 
minimize the residual variance. 



Contrary to the TLSE reconstructor, the MMSE reconstructor does not require any trun- 
cation. In fact, the MMSE filter includes prior knowledge of the SNR by the mean of the 



the numerator of Eq. 25 which weight the noise propagation. Thanks to this regularization 
term, the MMSE is able to optimize the reconstruction depending on the SNR: for frequencies 
that have a good SNR, MMSE and TLSE are equivalent, for poor signal to noise frequencies, 
the invert of the interaction matrix is weighted by the regularization term, avoiding noise 
amplification. It follows that, (i) even on a pure unseen frequency, the residual MMSE PSD 
is never higher than the signal variance, and (ii) whatever the truncation threshold used 
for the TLSE estimator, the MMSE approach always gives smaller residual variance (see 
also [4]). 

Following with the MMSE, it is interesting to note that for a frequency exactly on an unseen 
mode, the residual PSD writes as: 



Eq. 27 indicates that when the turbulence is equally distributed between the two layers 
(Ai=A2=0.5), the residual PSD exactly reaches the uncorrected spectrum. However, as soon 
as the turbulent strength profile is not uniformly distributed, the residual error is always 
lower than the uncorrected spectrum. We illustrate this behavior in Fig. 5 with three resid- 
ual PSDs computed for respectively (from top to bottom) [Ai=0.5,A2=0.5], [Ai=0.7,A2=0.3] 



regularization term a'^C^p^ (see Eq. 17). This regularization appears as an additive term in 




(27) 
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and [Ai=0.9,A2=0.1] (the [Ai=l,A2=0] case is superimposed with the on-axis reconstruc- 
tion error). This behavior can be explained by the fact that when the turbulence is equally 
distributed, the chance that the sum of both layers gives a null measurement is maximal. 
Statistically, no information can be extracted from the measurement. When the two layers 
do not have the same strength, this probability decreases and even for an unseen frequency, 
some signal can be measured. As a limiting configuration, if one layer is free of turbulence 
(Ai,2 = 0), we are able thanks to this prior information to entirely recover the incoming 
perturbation. As the MMSE filter includes the knowledge of turbulence strength per layer, 
it is able to discriminate these conditions, and to modulate the regularization term layer 
per layer. This also illustrates how the prior knowledge of the statistical conditions can 
consequently reduce the residual variance compared to the crude TLSE approach. 



0-^ 




Von Kormon PSD 

MMSE reconstruction [A, =0.5,A2 = 0.5J~ 
MMSE reconstruction [A, =0.8,Aj = 0.2] 
MMSE reconstruction [A, =0,9,A2 = 0. 1 ] 
On-Axis reconstruction 




Spotiol frequencies m 



Fig. 5. Influence of turbulence strength repartition on MMSE reconstruction 
error. 



5.C. Generalization to realistic cases 

To understand what is the impact of a larger number of layers/GSs on unseen frequencies, 
we first start with a generalization of the simple 2GS case. In Fig. 6 we show the residual 
PSD for the 2GS case, but treated as a 2-dimensional configuration. On the left, we keep 
the two layers atmosphere, whereas on the right we use the 10 layer profile as defined in 
Table 1. We only display here the results of the MMSE reconstructor. Note that a cut along 
the X coordinate on the left PSD of Fig. 6 would give the same MMSE reconstructor plot 
as the one presented in Fig.3. 
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Fig. 6. 2- dimensional residual PSDs for the 2GS constellation with the 2 layers 
profile (left) and the 10 layer profile (right). Log scale and inverted colors are 
used. The GS orientation is the same as in Table 2. 

If we compare the shape of the two PSD presented in Fig. 6, we easily understand what is 
the impact of multiplying the number of layers: unseen frequencies become sectors instead 
of localized energy peak. Indeed, with several layers, each pair of turbulent layer creates is 
own set of unseen frequencies, so the energy is spread upon more frequencies. 
The second important point to emphasize is that unseen frequencies only appear in the 
parallel direction of the GS geometry. As explained in Tokovinin & Viard [9], for all the 
frequencies perpendicular to the GS direction (/ = {fx=0;fy} for the 2GS case), the effective 
separation between the GSs in each layer (ai^2^n) will not change the measured phase of 
those frequencies. These frequencies are always perfectly measured, the noise propagation 
is the one of the classical AO situation. This explains the "clean" frequency area along the 
y coordinate in Fig. 6. The width of this "clean" frequency area only depends on the GS 
separation that dictates the lowest unseen frequency. 

We are now able to understand the generalization to the 4GS-8GS cases. In Fig. 7, we show 
the 2-dimensional residual PSDs for the MMSE (left), the TLSE (center), and the LSE 
(right) for respectively the 4GS (top) and 8GS (bottom) constellations. 

The 4GS case is nothing else than the product of two 2GS constellation rotated by 90°. 
We retrieve two "clean" frequency area oriented toward the GSs directions. The same 
argument can be applied to the 8GS constellation. In addition, the diameter of the 8GS 
constellation is larger than the 4GS one: the width of its "clean" frequency area is smaller. 

In Fig. 7, it also appears that the variance is much higher for LSE or TLSE reconstructors 
than for the optimal MMSE one as already explained for the 2GS case in Fig. 4. 

Finally, we note that as expected the number of low frequencies affected by unseen modes 
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Fig. 7. 2- dimensional residual PSDs for the 4GS (top) and 8GS (bottom) cases. 
From left to right: MMSE, TLSE, non truncated LSE. Log scale and inverted 
colors are used. All PSDs are displayed with same scale. The GSs orientation 
are the same orientation as in Table 2. 

is larger for the large FoV configuration that is the 8GS case than in the 4GS configuration. 
The noise level associated to each frequency is however lower because twice more GSs are 
used. It results that both the 4GS and the 8GS configuration give similar residual variance 
for the MMSE reconstruction: ~0.77 rd^ (with our simulation conditions). 

6. Model Errors 

Up to now, we have assumed that the priors were perfectly tuned, i.e., that they do corre- 
spond to real conditions. These priors comes from our knowledge of the system/atmospheric 
conditions, and are not always easily accessible. Moreover, they can evolve during an 
observation run and it could be impractical to re-compute the reconstruction matrix 
frequently. In any case, it is of prime importance to understand and quantify the impact of 
these statistical/model errors on the reconstruction process. 

Still using the three configurations presented in sect. 4, we first investigate the model errors 
with an error on the altitudes of the reconstructed layers and an error on the number of 
reconstructed layers. These errors occur both in the TLSE and the MMSE reconstructors. 
We then aim to (i) understand the impact of the model errors illustrated with the 2GS case; 
(ii) compare the relative performance and robustness of the TLSE and MMSE reconstructors 
with the 4GS configuration; (iii) compare the medium (4GS) and large FoV (8GS) systems 
in the MMSE approach. 
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6. A. Error on layer altitude 

We first start with the 1-dimensional 2GS case. We recall that only two layers are considered 
in this case. In Fig. 8 we show an example for which (i) the upper reconstructed layer is 
higher than real layer (Top panel) and (ii) the upper reconstructed layer is lower than real 
layer altitude (Bottom panel). In both situations, the residual PSD is strongly affected, and 
for some frequencies the residual variance rises above the uncorrected PSD. 
An error on the altitude where the layers are reconstructed affects the frequencies where 
unseen modes should be regularized/truncated. In fact, the altitudes of the reconstructed 
layers introduced in the model set the sensitivity of the reconstructors to unseen frequencies. 
A wrong geometry implies that wrong modes are going to be regularized/truncated. For 
the frequencies that correspond to well-seen frequencies for the model, the direct invert of 
the interaction matrix is performed, whereas these frequencies could correspond to badly 
seen modes in the real geometry. For theses frequencies, an over amplification of noise could 
appear. On the opposite, some frequencies that are well seen by the real geometry could be 
treated as unseen frequencies by the model. For these modes, whereas a direct inverse would 
have worked well, over regularization or truncation limits the accuracy of the reconstruction. 

In a more realistic configuration, the atmosphere includes several dominant layers, and 
this error will be minimized. Of course, the best match between real profile and dominant 
layers will give the best results. To illustrate this point, we investigate our realistic 10 layer 
cases (4GS-8GS). We reconstruct the 10 layers, but for each reconstructed layer we allow 
an error in altitude of ±X%, X going from 0% to 50%. For each altitude error, we perform 
50 random trials and we compute the mean residual variance. Results are shown in Fig. 9. 

We first compare the results obtained for the 4GS case. We find that both reconstructors 
(MMSE and TLSE) follow the same tendency: the residual variance increases with the error 
on the model altitude. In addition, the order of magnitude of this error is more or less the 
same for both reconstructors. 

If we now compare the results of the 4GS and 8GS cases, we notice that the large FoV 8GS 
configuration is much more sensitive to an error in layer altitudes than the medium FoV 
4GS one. For the medium FoV configuration, the impact of this error is moderated: an error 
of 50% on the layer altitudes increases the residual variance by ~20% with our turbulence 
profile. This indicates that a perfect knowledge of layer altitude is not necessary, as soon as 
many layers are reconstructed. However, this is no more true for the large FoV configuration: 
an error of 50% increases the residual variance by almost a factor of 2. As an analogy to 
what is explained by R. Ragazzoni [23], the larger the distance between reconstructed and 
true layer, the lower the maximal equivalent frequency reconstructed in the true layer. And 
this maximal equivalent cut-off frequency depends on the FoV. This can also be explained 
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Fig. 8. Residual PSD for an error in the altitude of the reconstructed layer. 
Top: Upper altitude is overestimated (8.5km instead of 8km). Bottom: Upper 
altitude is underestimated (7.5km instead of 8km). 



with Fig. 7: when the FoV is increased, the unseen frequencies cover a larger area and the 
noise amplification affects more and more frequencies. 

We conclude that this error term is impacting the tomographic reconstruction for large FoV 
configurations, and that either a good knowledge of the turbulence profile or the use of more 
GS to reduce the unseen frequencies area is required to limit the impact of this model error. 

6.B. Error on number of layers 

For computing reasons, the reconstruction is performed on a limited number of layers, gener- 
ally smaller than real turbulence conditions. In that configuration, the reconstructed volume 
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Fig. 9. 10 layers profile. Influence of an error on reconstructed altitudes of 
each layer. % of error means that each reconstructed layer is at an altitude of 
H±H*X%. 



does not match with real profile, and the resulting error can be important. We first illustrate 
the impact of an incorrect number of reconstructed layers with the 2GS case. As only 2 layers 
are used to model the real turbulence profile, a unique layer is used in the reconstruction. In 
Fig. 10, we illustrate this situation when the reconstructed layer is in the telescope pupil (full 
line) and when it is at an altitude of 4km (dashed-dotted line). It appears that the residual 
PSD is strongly influenced by this model error: the domain of frequency affected by noise 
amplification is significantly enlarged. Interestingly, when the reconstructed layer falls just 
between the two real layers, this noise amplification is reduced compared to a reconstruction 
in the telescope pupil. Indeed, we can show that if the real turbulence is equally distributed 
between the two layers, the optimal altitude for the reconstructed layer is just at the mid- 
dle. If the real turbulence is not equally distributed, the optimal altitude for reconstruction 
moves toward the strongest layer. 

We now investigate the effect of this model error in presence of a "real" turbulence profile. 
To do so, we use the 4GS/8GS case, the "real" turbulence profile is defined by the 10 
layers one, and we progressively increase the number of reconstructed layers. This has been 
introduced by Fusco [11] as Equivalent Layers (EL). We then follow their approach for the 
definition of the altitudes/strenght of the EL: the true profile is divided into N^i regularly 
spaced slabs, for each slab an equivalent height and strength is computed. Results are shown 
in Fig. 11. This figure shows that for the medium FoV constellation (4GS), a reconstruction 
on only few layers (typically 3 or 4) is enough to reduce consequently the remaining error. 
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Fig. 10. Residual PSD for a model error in the number of layers. The real 
profile is made of two-layers, whereas the reconstruction model use only one 
layer either in the telescope pupil (full line) or at an altitude of 4km (dashed 
dotted line). The MMSE and TLSE cases are superimposed. 

With the 4GS case, we also find that the MMSE residual variance converge more rapidly 
than the TLSE one. For the large FoV configuration (8GS), we observe that at least 8 layers 
must be reconstructed, and this number would certainly be larger if the initial profile had 
included more layers. As for an error on layer altitudes, the tomographic reconstruction is 
very sensitive to this model error, and particularly when the size of the FoV increases. We 
draw similar conclusions as in Sect. 6. A: if the number of GS is limited, a good knowledge 
of turbulence profile is necessary for an accurate reconstruction. 

7. Statistical Errors 

We have seen in the previous section the effect of the model errors, we now investigate the 
impact of the statistical errors. We then assume that the model is perfectly known: we always 
reconstruct the real number/altitudes of layers. Instead, we explore the consequences of (i) 
an error on the noise level (Cb 7^ Cb,t), (h) an error on the global strength of turbulence 
('^o 7^ ^o,t), and (iii) an error on the relative repartition of the turbulence (A„ 7^ A„,(). 

An error on these statistical parameters will impact on the value of the SNR used by 
the model (Cf. Eq. 17). Basically, the SNR used by the model sets the number of fil- 
tered/truncated frequencies. If this SNR is overestimated, the number of filtered/truncated 
frequencies decreases. For unseen frequencies, the model assumes that the SNR is good, a 
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Fig. 11. Influence of the number of layers used in the reconstruction process. 
Real profile is the 10 layers. 



direct invert of the interaction matrix is performed, and then we expect an over amplification 
of the noise localized on unseen frequencies. On the other hand, if the SNR is underestimated, 
the weight of the regularization term is increased: too much frequencies are filtered/truncated 
whereas they could have been corrected. We then expect a broadening of the PSD around 
the unseen frequencies, but no noise over-amplification. 

7. A. Error on noise priors 

In Fig. 12, we first use the 2GS case to illustrate the effect of an error on the noise estimation. 
The top plot in Fig. 12 shows the impact of an overestimation of the noise variance by 
a factor of 10 (Cb = 10Cb,t), and the bottom plot the impact of an underestimation of 
the noise variance by a factor of 10 (Cb = Cb,t/10). With this simple example, we indeed 
find that if the noise is overestimated (top figure), the residual PSD tends too rapidly to 
the uncorrected spectrum and too much frequencies are filtered/truncated. If the noise is 
underestimated (bottom figure), an over amplification of the residual error appears on the 
unseen frequencies and not enough frequencies are filtered/truncated. 

Although the residual PSDs are affected by a wrong noise model, at first glance, the impact 
of an error in noise statistics does not seem to increase significantly the residual error. To 
quantify this effect, we use our realistic cases (4GS-8GS) and we plot in Fig. 13 the residual 
variance for a MMSE/TLSE reconstructor having wrong noise priors. The error in noise 
priors is given in % of real noise: for instance, ±50% of error in noise variance corresponds 
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Fig. 12. Residual phase PSD for an error in noise priors. Top: Overestimation 
of noise by a factor of 10: Cb = lOCb.t- Bottom: Underestimation of noise a 
factor of 10: Cb = Cb,t/10. 



to Cb = ±1.5Cb,t- From Fig. 13, we learn that: 

• the residual variance in the MMSE approach is always smaller than the TLSE one, 
which means that even with wrong priors, the MMSE reconstructor gives better results 
than the TLSE one. 

• The TLSE reconstructor is slightly more sensitive to noise error than the MMSE one. 

• The large FoV configuration (8GS) is slightly more sensitive than the medium FoV one 
(4GS). 
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For both reconstructors, if the noise level is not perfectly known, a conservative ap- 
proach would be to overestimate the noise priors. 




Fig. 13. Residual variance of the MMSE(full-line) and TLSE (dashed-line) 
reconstructors with wrong noise model. Error in noise are given in % of real 
noise. Dotted-lines symbolize the minimal variance level when no errors on 
noise priors are committed. 



Another way to test the robustness of these reconstructors to the noise statistical error, 
is to study for a given reconstructor computed with a reference noise level, the impact on 
the performance when real conditions are changing. This is illustrated in Fig. 14 for the 4GS 
configuration. Both MMSE and TLSE reconstructors are computed for a given reference noise 
level ( "variation for real noise variance" equal 0), and we compute the residual variances when 
real noise conditions are changing. The change in real noise variance is given in % of noise 
priors: for instance, ±50% of variation in real noise variance corresponds to Cb,t = ±1.5Cb. 
In that situation, we find that the variation of the performance due to real conditions is much 
more impacting the residual variance than the noise model error itself. Indeed, an increase 
of the real noise variance by 50% leads to an increase of the residual variance by ~10%, 
whereas the corresponding model error only increases the residual variance by ~3% with our 
parameter set. Results for the 8GS configuration are very similar. 

Combining the results of Fig. 13 and Fig. 14, we conclude that a perfect knowledge of the 
noise priors is not required to obtain an accurate tomographic phase reconstruction. 
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Fig. 14. Robustness of the MMSE/TLSE reconstructor to changes of "real" 
conditions. The "real" noise level (defined by Cb,t) ranges from -90% to +90% 
around a reference level (0.5 rd^ when "% of deviation of the real noise vari- 
ance" = 0). The full and dashed lines show the residual variance for perfectly 
tuned MMSE and TLSE reconstructors: noise priors follow the real noise vari- 
ations. Dotted and dotted-dashed lines show the residual variance when noise 
priors are set to the reference noise level. 



7.B. Error on turbulent strength priors 
7.B.I. error on tq 

The global strength of the turbulence tq is included in the estimation of the SNR. The 
conclusions drawn in the previous section are then exactly transposable here: (i) the 
MMSE is always better than the TLSE with wrong truncation threshold (ii), the large FoV 
configuration is slightly more sensitive than the medium FoV one, and (iii) a conservative 
approach would be to underestimate the global strength of the turbulence. 



7.B.2. error on turbulence strength per layer 

The turbulence strength distribution is another essential parameter to be introduced in the 
regularization process. This term only affects the MMSE reconstructor, as the TLSE one 
does not include the profile distribution information. We have seen in Sect. 5.B that the 
turbulence strength distribution (A„) sets the SNR per layer. An error on the repartition of 
the turbulence strength will then produce an error on the SNR, and an over amplification of 
noise. 
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We illustrate this point in Fig. 15 (top) with a pessimistic example: the turbulence strength 
repartition used by the model is [Ai=0.1,A2=0.9], whereas real turbulence distribution is 
[Ai t=0.9, A2,t=0.1]. For this pessimistic situation, we see in Fig. 15 that the residual PSD 
rises above the uncorrected turbulent PSD. 




Von Karman DSP 
MMSE Reconstruction (tru. 
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Fig. 15. Error on turbulence strength priors with real turbulence distribution 
defined by [Ai,f=0.9, A2,t=0.1] . Top: turbulence strength repartition used by 
the model is [Ai=0.1,A2=0.9] (solid line). For comparison, we plot the residual 
PSD if the model were correct (dashed line). Bottom: turbulence strength 
repartition used by the model is [Ai=0.5,A2=0.5]. 



To avoid any noise amplification, a conservative approach would be to feed the model 
with a uniform turbulent strength distribution. This example is illustrated in Fig. 15 
(bottom). The true turbulent distribution is still [Ai^j=0.9, A2,t=0.1], the one used by the 
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model is [Ai=0.5,A2=0.5]. With such a model, we have seen with Eq. 27 that the residual 
PSD becomes independent of the modeled strength repartition. As expected, the constant 
profile is the one providing the less information. 

With the conditions of Fig. 15 (top), the total variance of the MMSE is higher than the 
TLSE one. However, following the conservative approach, i.e., the constant profile (Fig. 15 
(bottom)), the MMSE is always better than the best tuned TLSE. 

We could further investigate these model/statistical errors, for instance by combining 
several terms together. In that case, it is interesting to show that some model errors 
can be compensated from others (e.g. [24]). For instance, the global variance due to a 
mis-knowledge of layer altitude could be compensated by using more noise priors in the 
reconstruction. A pure quantitative study of model errors depends on system characteristics, 
and this exhaustive and specific work is certainly out of the scope of this paper. Instead, in 
the next section we choose to investigate the error due to the projection onto DMs. 

8. Projection onto DMs 

Up to now, we have only considered the tomographic phase reconstruction issue. In this 
section, we now investigate the projection onto deformable mirrors, and the supplementary 
error term related to generalized fitting. This term relies on Eq. 21 and depends on the 
optimization positions {{Pj}) and the number/altitudes of DMs. We mainly use the 4GS 
constellation, and we investigate different WFAO systems for a 42m telescope. We first 
consider a projection onto a single DM with one direction of optimization (MOAO/LTAO), 
we then increase the optimization field still assuming a single DM (GLAO), and finally we 
introduce several DMs with an optimization in the field (MCAO). 

Complementary to what have been done in the previous sections, we now extend the analysis 
to the whole field instead of concentrating on a particular direction. For each configuration, 
we choose to evaluate the performance in term of Strehl Ratio (SR), in 169 directions 
regularly spaced on a 13x13 grid across the field. Note that the SR metric is particularly 
relevant for diffraction-limited imagery (MOAO/LTAO/MCAO), but less appropriate for 
partial atmospheric corrections (GLAO, e.g. [25]). Nevertheless, for the sake of clarity and 
simplicity, we have decided to only consider SR in the following. 

For MOAO, optimization is done at the center of the field, for GLAO/MCAO optimization 
is done on the same directions. Fig. 16 illustrates the geometry used for the MOAO/LTAO 
configuration (top) and for the GLAO/MCAO configurations (bottom). Note that the GS 
constellation is rotated by 45° compared to Table 2. 

In all the following, except in Sect. 8. C. 3, no model/statitical errors are introduced. We 



31 



0.5 - 



-0.5 - 



-1,0 - 



o o 


actions 
direction 

O 


o 


o 


o 


o 


o 


o 


o 





o 





o 


o 


o 


o 


o 


o 





o 





o 





o 





o 


* 





o 


o 


o 





o 





♦ 





o 





o 








o 





o 





o 





o 





o 1 


o 


o 








o 


o 


o 





o 


o 


o 





o ; 


o 


o 








o 


o 


o 





o 


o 


o 





o 





o 


o 





o 





® 





o 





o 





o 





o 


o 





o 


o 


o 


o 


o 





o 





O - 





o 


o 


o 


o 


o 


o 





o 





o 





o : 


o 


o 


o 


o 


o 


o 


o 





o 


o 


o 





o 


o 


o 


♦ 


o 


o 


o 


o 


o 


o 


o 


* 


o 


o 





o 


o 





o 


o 


o 





o 





o 


o 


o 





o 








o 





o 





o 





o 





o r 



-1,5 I 

-1.5 -1.0 -0 5 0.5 -.0 1.5 



1,5 


□ Otpmisctio 
★ Guide Stor 
O PSF estim 


n di 


dire 
























1 ,0 




<5) 


© 


© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


© 






S) 


© 


© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


© 








<s) 


© 


* 


© 


© 


© 


® 


© 


® 


© 


* 


© 


© 




0,5 




@ 


® 


® 


© 


® 


® 


® 


© 


® 


© 


® 


© 


® 








© 




© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


® 








Si 


© 


© 


© 


© 


© 


© 


© 


© 


© 


© 


© 


© 




0,0 




© 


© 


© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


© 










© 


© 


© 


® 


© 


® 


© 


® 


© 


® 


© 


© 










® 


© 


© 


® 


® 


® 


® 


® 


© 


® 


© 


® 




0,5 




© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


© 


© 


© 








© 


© 


* 


© 


© 


© 


© 


© 


© 


© 


♦ 


© 


© 








© 


© 


© 


© 


© 


© 


® 


© 


© 


© 


© 


© 


© 




1,0 
1,5 




© 


© 


© 


© 


® 


© 


® 


© 


® 


© 


® 


© 


© 





-1.0 -0.5 0.0 0.5 '.O 

Aromin 



Fig. 16. Geometry used in the simulations. Top: One direction of optimization 
is considered at the center of the field. Bottom: optimization is performed all 
over the field. 

further assume that DM pitchs match WFS pitchs, which would correspond to DMs with 
~75x75 actuators. Atmospheric parameters are the ones introduced in Sect. 4. We stress 
that a different set of atmospheric parameters would lead to different numerical results. 
The derived performance is then only indicative and useful for a relative comparison of the 
systems. Finally, we recall that the telescope diameter is set to 42m. 

8. A. MOAO/LTAO: No projection error 

For some specific observations a uniform correction of the whole FoV is not required. For 
instance, in 3D spectroscopy extra-galactic studies, only few directions of interest must be 
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corrected for. This can be achieved by the way of MOAO or LTAO ( [26,27]). In these con- 
cepts, several GSs spread over the field are used to perform the tomographic reconstruction 
of the turbulent volume. The correction is then applied with one DM conjugated to the 
pupil per science channel: = Id. For LTAO, only one science channel is considered, for 

MOAO, few channels are corrected at the same time. For each direction of optimization j3, 
the projection onto the DM simply writes as: 

Popt = (28) 

In those directions, the performance is the one of the pure tomographic reconstruction 
described in Sect. 5, there is no supplementary term of error due to generalized fitting. An 
example of the expected performance of an LTAO/MOAO system is presented in Fig. 17 
(top) for the 4GS case. In this example, only one direction at the center of the field {(3=0) 
has been optimized. For this direction, the performance is only limited by unseen frequencies 
and it reaches ~50% of SR. Outside this direction, performance quickly decreases due to 
classical anisoplanatism. 



8.B. GLAO 

With still one DM conjugated to the telescope pupil, but an optimization in the whole 
field, we now investigate GLAO [28]. The GLAO correction represents the worst case for 
generalized fitting, as only one DM is used to correct the whole FoV. For a DM conjugated 
to the telescope pupil, P^^ simply becomes the identity matrix, and the projection term 
then writes as: 

Popt = {P%)p (29) 

it only represents an average over all the direction of optimization. An example of the 
expected performance with GLAO is presented in Fig. 17 (bottom). We see that the 
performance is limited to few % of SR, but it is almost uniform in the field, and the 
corrected area is larger than for the MOAO correction. For this example, the error due to 
generalized fitting at the center of the FoV can be as high as a factor of 5. 

8.C. MCAO 

To reduce the projection error, one must use more DMs optically conjugated to the layers 
in altitude. 



33 



30 




Fig. 17. SR maps for MOAO (top) and GLAO (bottom). GS are 45° rotated 
compared to Table 2. The SR have been hnearly interpolated on a 26x26 grid. 
For MOAO, the SR is ~50% at the center of the field. For GLAO, the SR is 
~10% at the center of the field. 



8.C.I. Generalized fitting 

When several DMs are used to perform the correction, the projection term is given by 

Eq.21. For a given optimized field, the generalized fitting error only depends on the number 

of DMs [12, 13]. The more DMs used to correct the volume, the better the match between 

mirrors position and turbulent layers, the better the performance in the field. 

To illustrate the effect of the generalized fitting error, in Fig. 18 we plot the mean SR over 

the field (computed over the 169 directions) as a function of the number of DMs. Error bars 

represent the standard deviation of the SR in the field. Altitudes of the DMs are defined by 

the altitudes of the equivalent layers, following the same procedure as in Sect. 6.B. Results 

are shown for the 4GS and 8GS cases using an MMSE reconstruction. 

Fig. 18 shows that with the 4GS (medium FoV) case, 3 DMs are sufficient to obtain almost 
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an optimal performance. However, for the 8GS (large FoV) case, this number increase to 8 
DMs. These results are very similar to those obtained in Fig. 6.B. Indeed, as introduced in 
Fusco et al. [4], the generalized fitting error is close to a model error on the number of layers. 




Fig. 18. Mean SR over the field in function of the number of DMs. Error bars 
represent the standard deviation of the SR. 

8.C.2. Impact of different reconstructors 

We now want to investigate the gain of the regularized reconstructor (MMSE) compared to 
the TLSE one. We use the 4GS case with 3DMs located respectively at [0, 3.5, 9.3]km to 
limit the impact of generalized fitting (see Fig. 18). In Fig. 19 we show the SR map for the 
TLSE (top) and MMSE (bottom) approach. We find the MMSE provides a performance ~1.2 
times better in mean SR: there is a significant gain to use a regularized algorithm compare 
to a crude TLSE one. Particularly in the corners of the field where the MMSE extrapolate 
the phase estimation, whereas the TLSE becomes very sensitive to unseen frequencies. Note 
that in the TLSE case, another choice of threshold could make the performance better on 
the GS (as good as the MMSE reconstruction), at the price of decreasing the performance 
in the field. 

8.C.3. Equivalent Layers 

Another classical approach in MCAO is to reconstruct the turbulent volume only at the 
DMs altitudes. No projection onto the DMs is then required and Popt = Id. Following our 
approach, this is nothing else that a model error on the number of layers as described in 
Sect. 6.B. In Sect. 6.B we focused on the effect of this error at the center of the field, we 
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Fig. 19. SR maps for an MCAO working with 3DMs at [0, 3.5, 9.3]km. Com- 
parison between TLSE (top) and MMSE (bottom) reconstructors. SR[min, 
max, mean, standard deviation] = [10%, 42%, 25%, 6%] for the TLSE and 
[30%, 57%, 47%, 6%] for MMSE 



now want to evaluate the impact of this approach for the performance in the whole field. We 
choose to simulate a system working with 3DMs as described in Sect. 8. C. 2. In Fig. 20 we 
compare the respective performance obtained for the TLSE (top) and the MMSE (bottom) 
reconstructors. 

Similarly to the results of Sect. 8.C.2, we first find that the MMSE approach provides 
a performance ~1.2 times better than the TLSE one. Another interesting exercise is to 
compare the results of Fig. 20 with those of Fig. 19. Doing so, we first find that the mean 
SR are very close for both approaches: a reconstruction performed directly on the DMs only 
gives a 5% less performance in SR. This is indeed consistent with the results obtained in 
Fig. 11 for the center of the field. Note that this errors strongly depends on the turbulence 
profile, and different turbulent conditions would give different results. In any cases, the main 
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Fig. 20. SR MMSE in the equivalent layers approach: turbulent volume is 
reconstructed only at DMs altitudes. Top TLSE reconstructor SR[min, max, 
mean, standard deviation] = [7%, 47%, 21%, 7%] Bottom: MMSE reconstruc- 
tor. SR[min, max, mean, standard deviation] = [24%, 58%, 42%, 6%] 



difference between Fig. 20 and Fig. 19 is that the correction is no more uniform over the 
field when one uses a reconstruction directly onto the DMs. Indeed, the optimization of 
the performance for specific directions comes from the projection term Popt- To enable an 
optimization of the performance in the field, one must then reconstruct the turbulent volume 
on more layers than DMs. The other main draw-back of a reconstruction directly on the DMs 
is that the corresponding model error is very sensitive for large FoV systems (see Sect. 6.B). 
For instance, the same study for the 8GS configuration only gives a mean SR of 15% for the 
MMSE reconstructor. The main advantage of the EL approach is that the tomography is 
much more simplified compared to the full MCAO presented in Sect. B.C. 2. 
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9. Conclusion 



We have first presented a matrix formalism for the Fourier modeling of any WFAO systems. 
Our Fourier tool includes all the specifies of the WFAO systems such as tomography, num- 
ber/positions of DMs, model/statistical errors. The Fourier approach is interesting because 
it allows a fast and easy exploration of the broad parameter space, as well as a detailed 
comprehension of the underlying physical phenomena. For ELTs studies, it offers a fine and 
accurate modeling tool, able to provide the end-product PSFs. 

Based on this Fourier tool, we have explored three main issues of any WFAO system, respec- 
tively: unseen frequencies, model/statistical errors and projection errors. Our goal was not to 
derive numerical results, but rather to point out general trends shared by the tomographic 
systems. We first illustrated how the GS geometry set the amount of unseen frequencies. 
The repartition of unseen frequencies is essential to understand the additional errors: they 
represent the roots where additional errors are growing. 

Our exploration of model/statistical errors draw several limitations shared by all WFAO sys- 
tems. First, we retrieve that the MMSE approach is always more accurate and robust than 
the best TLSE one. Then, that the tomographic reconstruction is robust to statistical errors, 
but is very sensitive to model errors. For medium FoV systems, a good knowledge of any 
statistical/atmospheric priors is not essential. However, for large FoV systems the sensitivity 
to turbulent distribution errors becomes significant. This implies that large FoV concepts 
would require a large number of GS to limit unseen frequencies, or a good monitoring of 
turbulent conditions. 

We also derive few rules for a robust control when system and atmospheric conditions are 
not perfectly known. In such cases, it is advice to (i) overestimate noise variance, (ii) under- 
erestimate turbulence strength, (iii) tends to a constant turbulence profile strength, (iv) use 
more layers in the reconstruction. Finally, we give a first insight into typical performance of 
the future WFAO systems for ELTs, and we show that regularized tomographic algorithms 
are essential to provide a significant gain in performance for these future systems. 
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